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We compare two recently proposed methods for the characterization of phase transitions in 
small systems. The usefulness of these techniques is evaluated for the case of structural transition 
■ in alanine-based peptides. 
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1. INTRODUCTION 

When phase transitions are studied in statistical physics, the basic assumption is usually that the dimensions of 



, the macroscopic system are very large compared with that of the constituting elements. However, there are many 
00 ' 

phenomena in finite systems which resemble phase transitions, for instance, in the physics of clusters of atoms 



The main question is how the observed effects in small systems can be related to true phase transitions in 



macroscopic (or infinite) systems. Attempts in this direction include the exploration of the density of complex zeros 

^ ' ..... a 

of the canonical partition function for finite and small systems by Borrmann ct. al. |2j] and the linear behavior for 

this limiting density HQ . The approach suggested by Janke and Kenna ( JK) presents also a new scaling relation 
Eh ' to identify the order and strength of a transition from the behavior of small systems ||. 

In this paper we try to evaluate the usefulness and validity of the Borrmann and the JK approaches for the 

£ '■ 

investigation of structural transitions in biomolecules, another important example for "phase transitions" in small 

O ■ 

systems. For this purpose we consider the helix-coil transition in polyalanine and test whether the two approaches 
allow a characterization of that transition. In addition, we also investigate a second molecule Alaio-Glys-Alaio in 



the context of Borrmann approach. 



2. ANALYSIS OF PARTITION FUNCTION ZEROS IN SMALL SYSTEMS 



In the canonical ensemble a system is completely described by its partition function 
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Z((3) = V n (£)exp(-^). (1) 
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Introducing the variable u = exp(— kf3) with conveniently defined constant k allows to write the partition function 
for discrete energy models as a polynominal: 

Z(u)=J2ME)u E . (2) 

E 

We label the complex zeros in order of increasing imaginary part as Uj = exjp(—k(3j),j = 1,2, .... 

As established by Yang and Lee || and later by Fisher 0, the statistical theory of phase transitions can be 
described by the distribution of complex partition function zeros. In the case of a temperature driven phase 
transition, we expect the zeros (or at least the ones close to the real axis) condense for large enough system size L 
on a single line 

Uj = u c + rj exp(i(p) . (3) 

As the system size L increases, these zeros will move towards the positive real w-axis and for large L the corre- 
sponding value is j3 c , the inverse of the physical critical temperature T c . Crucial information on phase transitions 
can be obtained from the way in which the first zero approaches the real it-axis. However, such an analysis depends 
on the extrapolation towards the infinite large system and does not allow characterization of the situation in small 
systems. 

One possible extension of the above ideas to "phase transitions" in small systems is the classification scheme by 
Borrmann et al. Q . This classification scheme computes density of zeros taking into account numerical estimates 
of the first four complex zeros. Now, writing the complex zeros as u k = Re(ufc) + i Tk, where Tk stands for Im(itfc), 
the assumed distribution of zeros on a straight line allows to define two parameters a u and 7 U : 

ln^(r 3 )-ln^) 1 / 1 1 \ 

a-u = : : , where 0(t*) = - ] - + , r , (4) 

lnr 3 -lnr 2 2 \\u k - u k -i\ \u k+ i ~ u k \ ) 

with k labelling the first zeros, and 

7„ = (Re(« 2 ) - Re(ui))/(r 2 - n) . (5) 

Note, that our notation differs from that in Ref. || in that we define the discrete line density <p m function of the 
it-zeros instead of the inverse temperature (3. Following the classification scheme by Grossmann and Rosenhauer 
phase transitions can now be classified according to the values of these two parameters: for a u — and j u = 
one has a phase transition of first order, it is of second order if < a u < 1 and arbitrary 7„, and for a u > 1 and 
arbitrary "f u one has a higher order transition. Of course, as expected for a phase transition, the imaginary part of 
the leading zero (ti) should move towards the real axis to obtain in the thermodynamic limit a real temperature. 
However, for small systems with finite ti, we may have a < 0, which translates the condition of first order 
transition to a u < and 7„ = 0. This classification of phase transitions has been tested for the finite Bose-Einstein 
condensates in a harmonic trap, small magnetic clusters and nuclear multifragmentation 

Another extension of partition function zero analysis to small systems is the approach by Janke and Kenna (JK) 
|^| , which uses that the theoretical average cumulative density of zeros 

G dr 3 ) = ^ (6) 
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can be written in the thermodynamic limit and for a first order transition as 

Goo(r) = 9oo (Q)r + ar w+1 + ■ ■ ■ . (7) 

Here, the slope at the origin is related to the latent heat, Ae oc goo(0). Equations (||) and (0) imply that the distance 
r\j of a zero from its critical point can be written for large enough lattice sizes as Im Uj{L) since Tleuj(L) ~ u c . 
Hence, in this limit Eq. (||) and (Q) lead to the following scaling relation for the cumulative density of zeros as an 
equation in j and L, 

2j-l 



2L d 



«i(Im Uj (L))" 2 +a 3 . (8) 



A necessary condition for the existence of a phase transition is that 03 is compatible with zero, else it would 
indicate that the system is in a well-defined phase. The values of the constants a\ and 02 then characterize the 
phase transition. For instance, for first order transitions should the constant 02 take values 02 ~ 1 for small r, and 
in this case the slope of this equation is related to the latent heat through the relation 0] 

Ae = ku c 2ttcli , (9) 

with u c = exp(— kp c ). On the other hand, a value of ai larger than 1 indicates a second order transition whose 
specific heat exponent is given by a = 2 — 02. 

The above approach was originally developed and tested for systems with well defined first order phase transitions 
such as the 2D 10-state Potts and 3D 3-state Potts model. The obtained results agree with previous work from 
numerical simulations and partition function zeros analysis of systems up to L = 64 |^|| (L = 36 for the 3D case 
]To|]). However, some difficulties in the identification of the order of the transition appear when applied to a model 
with a very weak first order transition [ pT[ . 

Here we explore further these two approaches for a non-trivial example of structural transitions in biomolecules, 
namely the helix-coil transition in polyalaninc. The characteristics of this so-called helix-coil transition have been 
studied extensively fL2|| . For instance, evidence was presented |]l3|-|l5| that polyalanine exhibits a phase transition 
between the ordered helical state and the disordered random-coil state when interactions between all atoms in the 
molecule are taken into account. These previous results, obtained independently by different methods, are used 
here as benchmarks for the test of the usefulness and to point out the difficulties in applying the classification 
scheme of Borrmann and Janke/Kenna for research of structural transitions in biomolecules. 

Our investigation of the helix-coil transition for polyalanine is based on a detailed, all-atom representation of that 
homopolymer. Since one can avoid the complications of electrostatic and hydrogen-bond interactions of side chains 
with the solvent for alanine (a nonpolar amino acid), explicit solvent molecules were neglected. The interaction 
between the atoms was described by a standard force field, ECEPP/2, Chains of up to N = 30 monomers 

were considered, and our results rely on multicanonical simulations Jl7| of N sw Monte Carlo sweeps starting from 
a random initial conformation, i.e. without introducing any bias. Our statistics consists of iV SUJ =400,000, 500,000, 
1 ,000,000, and 3,000,000 sweeps for N = 10, 15, 20, and 30, respectively. Measurements were taken every fourth 
Monte Carlo sweep. Additional 40,000 sweeps (N = 10) to 500,000 sweeps (N — 30) were needed for the weight 
factor calculations by the iterative procedure described first in Ref. jlffl. In contrast to our first calculation of 
complex zeros presented in Ref. |l4|| , where we divided the energy range into intervals of lengths 0.5 kcal/mol in 
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order to make Eq. (|^) a polynomial in the variable it = e - ' 3 / 2 , we avoided any approximation scheme in the present 
work. This because the above approximation works very well for the first zero, but not for the next ones. Since 
we need high precision estimates also for the next zeros we applied the scan method (see Ref. [Q and references 
therein) . 

In Table 1 we present our first four partition function zeros, although the fourth one is less reliable due to the 
presence of fake zeros Jl9[ ]. Using these zeros we first calculated the parameters a u and 7„ which characterize in 
the Borrmann approach phase transitions in small systems. As one can see in Table 2 the so obtained values for 
polyalanine are characterized by large fluctuations. It seems that the median of the a u values is a u = which 
would indicate a first order transition. However, our data are not good enough to draw such a conclusion. Our 
error estimates were obtained by means of the Bootstrap method J2(^,[ll| . For this reason, we tried instead the JK 
scaling relations. Table 3 lists the parameter a3(N) of Eq. (|s[> . Here, the average cumulative density of zeros is 
replaced by 

Gn(tj) = ^^, (10) 

where we have translated the linear length L as N 1 ^ fL4| . Therefore all finite size scaling relations can be written 
in terms of the number of monomers TV. In Fig. 1 we show the cumulative distribution of zeros in a log-log plot 
for chain lengths N = 10, 15, 20 and 30. 

The values 03 are compatible with zero for chains of all length indicating that we have indeed a phase transition. In 
order to evaluate the kind of transition we also calculate the parameters a± (N) and a 2 (N) which we also summarize 
in Table 3. Our analysis procedure differs from JK in that we are calculating a 2 for a fixed iVbefore performing a 
finite-size scaling analysis of that quantity, i.e. we account for the dependence of the average cumulative density on 
the system size. The parameter a 2 (N) decreases with system size and a log- log plot of this quantity as a function 
of chain length, as shown in Fig. 2, suggests a scaling relation 

a 2 {N) = a 2 + bN- c . (11) 

A numerical fit of our data to this function leads to a value of a 2 = 1.31(4) with goodness of fit Q = 0.95. Using 
a = 2 — a,2 we find a = 0.70(4) which is barely compatible with our previous value of a — 0.86(10) in Ref. [fl4|| , 
obtained from the maximum of specific heat. A fit of all four chain lengths can also not exclude a value a 2 = 1 
(i.e. a first order phase transition) since we can find acceptable fits with Q > 0.55 in the range 0.92 < a 2 < 1.44. 
Especially, a close examination of Fig. 1 and Fig. 2 shows that the N = 30 data point exhibits a considerable 
deviation from the trend suggested by the smaller chain lengths. Since the N — 30 data are the least reliable, we 
also evaluated Eq. jll| ) omitting the N = 30 chain. This leads to a value of a 2 — 1.16(1) and a critical exponent 
a = 0.84(1) which is now compatible with our previous result a — 0.86(10). 

While the JK approach is able to reproduce results for polyalanine obtained in previous work fl4[] , it does not 
allow to establish the order of the helix-coil transition from simulation of small chains. Our results for the parameter 
a 2 seem to favor a second order transition, but are hampered by large errors (such that a first order transition 



cannot be excluded) and disputed in Ref. |15| were indications for a finite latent heat were found. 

Borrmann's approach led to even less decisive results for polyalanine since the fluctuation were too large to draw 
firm conclusions. However, the advantage of this approach is that it does not require a finite-size scaling. We have 
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therefore used this method to investigate a second molecule, the slightly more complicated Alaio-Glys-Alaio. A 
detailed study of this molecule will be published elsewhere |2l]] . 

We describe the interactions between the atoms again by the ECEPP/2 (ll| force field (as implemented in the 
program package SMMP |2^]). Our results rely on multicanonical simulations of 4, 000, 000 Monte Carlo sweeps 
starting from a random initial conformation, i.e. without introducing any bias. Measurements were taken every 
tenth Monte Carlo sweep. Additional 500,000 sweeps were needed for the calculations of the weight factor. 

Analyzing the partition function zeros for this peptide, we find a "phase transition" at two temperatures, each 
being characterized by a line of complex zeros. The corresponding first four zeros for each characteristic line 
are listed in Table 4. Both temperatures also correspond with peaks in the specific heat which is displayed in 
Fig. 3 as a function of temperature. In Ref. |2l[] it is shown that the higher temperature T\ — 480 K marks the 
helix-coil transition, i.e. the temperature where secondary structure elements (a-helices) are formed. At the lower 
temperature Ti = 260 K, the peptide then folds into its native structure, a U-turn-like bundle of two (antiparallel) 
alpha-helices connected by a turn. Here, we use the Borrmann approach to study these two transitions in more 
detail. Using the zeros from Table 4 we calculated the parameters a u and j u which characterize in the Borrmann 
approach phase transitions in small systems. For the first transition, at T = 480 K, which corresponds to the 
above studied helix-coil transition in polyalanine, we find a u = —1.46 and j u = 0.49. While the negative value of 
a u indicates a first order transition, our results do again not allow us to draw a firm conclusion on the nature of 
the transition since the values of the two parameters a u and "f u are subject to large fluctuations when we vary the 
number of sweeps considered in the calculation of the zeros. This is not surprising since we were also not able to 
establish clearly the order of the helix-coil transition in in the previous example of polyalanine. However, our data 
are more decisive in the case of the second transition, at T — 260 K, which marks the compactification and folding 
of the peptide. Here we find a u = 0.38 and 7 U = —1.07. These values indicate a second-order transition which is 
consistent with what one would expect for a transition between extended and compact structures. We hope that 
by doubling the number of Monte Carlo sweeps we will be able to verify this results and also obtain reliable values 
for the Borrmann parameter characterizing the helix-coil transition at T = 480 K. These simulations are now under 
way. 

3. CONCLUSION 

Let us summarize our results. We have evaluated two recently proposed schemes for characterizing phase tran- 
sitions in small systems. Simulating polyalanine molecules of chain lengths up N — 30 residues by multicanonical 
Monte Carlo, we calculated the partition function zeros of these molecules. Analyzing these zeros by the JK 
approach we were able to reproduce for polyalanine results obtained in previous work ]l4| while the Borrmann 
approach led here to inconclusive results. However, the later method seems more appropriated to study structural 
transitions in proteins since it does not require a finite size scaling. When applied to Alaio-Glys-Alaio, the Bor- 
rmann approach allowed us to characterize structural transitions in this molecule. However, our results are limited 
again by the precision with which one can calculate the first four partition function zeros. While these numerical 
limitations restrict the applicability of the Borrmann and the JK approaches, our results demonstrate that both 
methods can be useful tools for the study of "phase transitions" in biomolecules. 
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Figure Captions: 



Figure 1. Distribution of zeros for the polyalanine model. 

Figure 2. The parameter 02 (iV) for polyalanine molecules of length N in a log- log plot. 
Figure 3. The specific heat C(T) as a function of temperature for Alaio-Glys-Alaio- 

Table Captions: 

Table 1. Partition function zeros for polyalanine. 

Table 2. Estimates for the parameters a u and j u for polyalanine. 

Table 3. Estimates for the JK parameters for polyalanine. 

Table 4. Partition function zeros corresponding to the two transitions observed for Alaio-Glys-Alaio- 
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Table 1. Partition function zeros for polyalanine. 


N 


Re(ui) 


Im(ui) Re(u 2 ) Im(u 2 ) 


Re(u3) Im(w 3 ) 


Re(u 4 ) 


Im(u4) 


10 
15 
20 
30 


0.30530(12) 
0.356863(61) 
0.374016(41) 
0.378189(19) 


0.07720(14) 0.2823(13) 0.13820(61) 
0.053346(39) 0.34167(60) 0.10440(59) 
0.042331(45) 0.36161(27) 0.08109(24) 
0.027167(32) 0.37399(14) 0.05420(27) 


0.2459(72) 0.1851(63) 
0.3331(48) 0.1454(28) 
0.3569(27) 0.1154(13) 
0.3693(11) 0.0804(13) 


0.172(11) 
0.3067(81) 
0.3336(56) 
0.35854(63) 


0.2200(71) 
0.1689(32) 
0.1470(27) 
0.1022(43) 






Table 2. Estimates for the parameters 


a u and 7„ for polyalanine. 






N 




a« 








10 
15 
20 
30 




-0.36(17) 
0.41(19) 
0.06(14) 
0.19(14) 






-0.365(17) 
-0.291(11) 
-0.3229(78) 
-0.1568(58) 


Table 3. Estimates for the JK parameters for polyalanine. 


N 




ai 








10 
15 
20 
30 




6.17(60) 
4.37(19) 
3.62(26) 
3.54(31) 


1.862(46) 
1.664(16) 
1.558(24) 
1.473(30) 




0.01(14) 
0.014(69) 
-0.014(98) 
-0.007(61) 
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Table 4. Partition function zeros corresponding to the two transitions observed for Alaio-Glys-Alaio. 



Re(tti) 


Im(tti) 


Re(u2) 


Im(ti2) 


Re(u 3 ) 


Im(w 3 ) 


Re(u4) 


Im(u4) 


0.349620 


0.047937 


0.361590 


0.072467 


0.320531 


0.097282 


0.370090 


0.132103 


0.144820 


0.039763 


0.104362 


0.077531 


0.074312 


0.100345 


0.033091 


0.116098 
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Fig. 1. Distribution of zeros for the polyalanine model. 
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Fig. 2. The parameter a,2{N) for polyalanine molecules of length N in a log-log plot. 
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Fig. 3. The specific heat C(T) as a function of temperature for Alaio-Glys-Alaio. 
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